ANALYSIS OF EMISSION CONTROL BY A CATALYTIC RECUPERATOR AND BY 

UPSTREAM MONOLITHS 



INTRODUCTION 




sing post turbine 

catalytic monoliths and catalytically active exhaust-side recuperator surfaces to remove 
hydrocarbon and carbon monoxide emissions from the Capstone turbogenerators- 




lis report represents the detailed analysis, 
calculations, results, and conclusions. Also attached are the MathSoft Mathcad*^^ 6+ documents 
comprising the analyse;s arid an Excel 5.0 spreadsheet of the findings. 



TRANSPORT AND KINETIC ANALYSIS OF THE RECUPERATOR AND 
MONOLITHS 

The analysis of the conversion of heptene (C7H14), methane (CPLt), and carbon monoxide 
(CO), were completed for several configurations of a recuperator and honeycomb catalytic 
converters. The method used was a MATHCAD program, and the procedures were similar to the 
document sent to Capstone^^^|||||||||||^CAFSTONE.MCD*document)^ 

The analyses consist of the following features (see printouts-for-GPSTQNE^5^G D and 
-ePST0NE6.MGD): 

Gas Properties 

• The properties of the gases used in the combustion mixture are described on 
pages 1-8. 

• The gas compositions both for complete combustion air with fiiel equivalence 
ratio 0.1 in moist (2 vol% water vapor) air, as specified in the problem 
statement, and also for an adiabatic temperature rise fi-om 1 000° to 1 5 1 1_F 
using heptene are defined on page 2. The diflference in properties between 
these two compositions is insignificant. 

• Essential transport properties, such as the diffusivities of the added gases are 
determined on pages 3-4. 

• Mixture properties such as viscosity, thermal conductivity, heat capacity, and 
various dimensionless transport numbers such as the Prandtl and Schmidt 
numbers are defined and calculated on pages 5-8. 



Recuperator and Converter Dimensions and Operating Conditions 

• The dimensions of the honeycomb sizes from 100 to 800 cells/in.2 are defined 
on page 9 along with the size and dimensions of the recuperator and the 
upstream monolith converters. 

• The five operating combustion conditions representing the different power 
level cases for the turbogenerator are specified on pages 9-10. 

• The pressure drop for the 600 or 800 cells/in^ recuperator is determined on 
page 1 1, and the analysis of heat transfer within the recuperator is described on 
page 12. 

Recuperator Performance 

• Mass transfer rates to the wall of the recuperator (assuming very fast 
conversion kinetics at the wall) are described on pages 13-14 for all three 
added components (C7H14, CH4, and CO). These calculations used average 
temperatures to determine an integrated formula for fuel conversion. 

• An alternative form of conversion was developed by looking at local Nusselt 
number (modified for the square channel configuration). This analysis is 
described on pages 15 and 16. The two methods show good agreement for 
C7Hi4Conversion. 

• Catalytic properties and fiinctions are defined on pages 17-18, and the specific 
rate expression for oxidation of dilute CH4 by supported palladium oxide 
(PdO) is described on page 19. 

• Slow (in CPST0NE5.MCD and related to the rate of CH4 oxidation) and 
moderate (in CPST0NE6.MCD and related to the rate of C3H6 oxidation) rate 
constants for C7H14 oxidation are described on p. 19. 

• The performance of the recuperator for conversion with temperature-dependent 
mass transfer properties and specific catalytic rates for C7H14 and CH4 are 
described on pages 20-2 1 . 

Performance of Upstream Monoliths 

• The operating conditions for 200, 400, and 600 cells/in.^ honeycombs is 
described on pages 22-25. 

• The pressure drop for three honeycombs (200, 400, and 600 cells/in.^) is 
described on pages 22-23, and matched to the maximum specified conditions. 
From of these calculations, two honeycomb converter configurations, 400 
cells/in.2 by 4 in. length and 600 cells/in.2 by 3.5 in. length, were chosen for 
mass transfer calculations. 

• The effluent concentrations of added components are calculated for the 
constant temperature monoliths (400 cells/in.2 in CPST0NE5.MCD and 600 
cells/in.2 in CPST0NE6.MCD) for the five operating conditions (cases) on 
pages 24-25, assuming very fast catalytic rates. 

• Finally, on pages 25-27, the conversion of C7H14 and CH4 using finite catalytic 
rates are described for the honeycomb converters (400 cells/in.^ in 
CPST0NE5.MCD and 600 cells/in.2 in CPST0NE6.MCD). 




RESULTS OF THE ANALYSIS 

The principal results of the analysis of post turbine fuel component conversion are shown 
in the attached spreadsheet tables (CAPSTN5.XLS) and Mathcad documents (CPST0NE5.MCD 
and CPST0NE6.MCD). The main results of the analyses are summarized as follows: 

Catalytic Recuperator 

• The AP/P calculations seem to ve ry closely fit your measured values. 4HBH^ 
ft/ll/ltKtf/f^ We chose as the best calculation of pressure 

oropai^OO cells/in^xT-m^ith 10% additional added length to describe the 
chevron geometry of the recuperator foils. 

• The estimated averaged wall temperatures, or ski n temperatures, within the 
recuperator seemed a little low (< SO^F) based on^fl^BHHlHBHHI^ 

discussions of recuperator thermal efficiency. As a consequence, the cell size 
was increased to dOO cells/in^ keeping the length the same (7 in. x I.l) for a 
conservative calculation of mass transfer rates. 

• For the mass transfer limited cases (sissuming fast kinetics), virtually complete 
conversion CH4 and CO is predicted (see Table of Results 1). Except at the 
highest power level, conversion of CH4 and CO was essentially complete 
within the first inch or two of entry into the recuperator. 

• For the transport-limited cases with C7H14, maximum conversions are low 
only for Case 5 (full power). Even in this case, 50 ppm would be converted to 
< 1 ppm (0.2 ppm) while concentrations entering the recuperator above 250 
ppm would exit at levels above 1 ppm 

(500 ppm => 2,1 ppm). The situation improves for the other cases (e.g.. Case 4 
gives 500 ppm => 0.02 ppms while 5000 ppm => 0.17 ppm) . 

• Conversions would be much higher for transport limited conversion in a the 
800 cells/in.2 recuperator, as calculations with this cell size were found 
reproduce the observed pressure drop. In this case, the logarithm of the inlet 
concentration divided by the outlet concentration for all added components 
would be '-'35 % greater with the smaller cell size. 

• The step-wise calculation which includes the temperature on flow velocity and 
transport properties showed reasonably good agreement with the integrated 
calculation xising average properties. This method gave lower overall 
conversion than the averaged property formulas, but the results were close. 
Perhaps the higher property values at the inlet, particularly the greater velocity, 
more than offset the lower transport values at the outlet side. 

• CH4 conversion is catalytically rate limited under all conditions (see Table of 
Results 1). Conversions were calculated to be very low, with a maximum 
conversion of -66% under Case 1 and -32% conversion under Case 5. These 
conversions would increase somewhat, perhaps by a factor of 2, with the use of 
a thicker coating (e.g., 25 fim vs. 10 |im). 

• CO conversion is likely to exceed 99% in all operating conditions for the 
recuperator. This is because the kinetics of CO conversion are probably much 
faster than CH4 and may well be close to the transport-limited case, at least 
through the (hotter) first several inches of the recuperator where most of the 
CO would be consumed. 



• The catalytic recuperator, C7H14 rate determines its conversion in the recuperator 
conditions (see Table of Results 1). For the slower rates (derived from the rates 
for CH4 oxidation) using the 600 cells/in ^ value for recuperator cell size, we 
find a range from 99.4% conversion under Case 1 to about 85.6% under Case 5. 
If twice the rate of combustion of propane is extrapolated to the higher 
temperatures of the recuperator, much higher C7H14 conversions result, 99.98% 
under Case 1, and 96.4% under Case 5. 

Catalytic Monoliths 

• The pressure drop reaches maximum specifications for the 400 cells/in.^ 
honeycomb at 4-in. length. Maximum pressure drop would be reached in a 
shorter (-3.5 in. length) 600 cells/in.^ monolith. Both configurations were 
chosen for analysis (CPSTONE5.MGD for 400 cells/in.^ and 
ePST0NE6.MGD for 600 cells/in 2). 

• Transport-limited conversion shows excellent performance for methane and 
CO in all cases, including the fiill-power case, for both honeycomb 
configurations (see Table of Results 2). 

• As for the recuperator, methane conversion is again reaction rate limited (see 
Table of Results 3), with conversions range fi'om 93.8% under Case 1 with the 
600 cells/in.2 monolith to only 31% conversion under Case 5 for the 400 
cells/in.2 monolith. 

• C7H14 conversions under transport-limited conditions depends on the specific 
configiu-ation (see Table of Results 2). For the 600 cells/in.2 monolith 
conversions are very high for Case 1 , but fall to about 92% for Case 5. With 
the 400 cells/in.2 configuration, these fast kinetic transport-limited conversions 
fall to 87% conversion under Case 5. 

• The reaction rate of C7H14 strongly affects its conversion in either honeycomb 
monolith (see Table of Results 3). For the extrapolated propane kinetics at 600 
cells/in.2, 98% conversion is expected under Case 5, whereas the 400 cells/in.2 
using the more conservative extrapolation using PdO methane kinetics, 86% 
conversion is expected under Case 5. 

Effect of Washcoat Thermal Conductivity in a Recuperator 

• The application of a thin washcoat of porous al umina gives a v ery slight 
increase in the overall heat transfer coefficient^M|||MH^ 

ePSNALT2.MCD). The thermal conductivity of the washcoat, while 
significantly less than that of stainless 321, is greater than that of the gas phase 
boimdary. Because the washcoat displaces the gas phase (creating a narrower 
channel) it causes a slight increase in overall heat transfer coefficient. 



OVERALL CONCLUSIONS AND RECOMMENDATIONS 



The conclusions and recommendations drawn from the resuhs of the analysis is listed 

below: 

1) The actual C7H14 outlet concentration from the turbine are needed to assess the 
effluent emissions in the high-power case (Case 5) for either control strategy. 
Provided that fast catalytic rates can be provided in washcoats on either the 
recuperator or upstream monoliths, maximum conversions for C7H14 under fiill 
power range from 99.6% (for an effective cell size of 600 cells/in^) to > 99.9% 
(if the effective cell size is 800 cells/in^). 

2) The two methods (averaged properties and integrated formulas vs. local 
temperature dependent properties and finite difference integration) show good 
agreement for recuperator performance with infinitely fast catalytic reactions at 
the wall. 

3) Clearly, it will be necessary to apply and measure specific catalytic rates for 
C7H14 or gasoline to be sure that these estimated conversions could be 
achieved under actual operating conditions for either the catalytic recuperator 
or the upstream honeycomb monoliths. 

Thin wash coats (< 20 ^m) should have minimal effect on the transport of heat 
through the recuperator given the much greater resistance of the gas phase 
boundary layers. 

For the current operating conditions, use of a standard 400 cells/in^ x 
4-in length catalytic monolith upstream of the recuperator is recommended. 

6) If the recuperator temperature can be raised 200*"?) high conversions of CO, 
C7H14, and even CH4 can be expected under all flow conditions. In this case it 
is possible that only the first 3 to 4 inches of the recuperator be wash coated 
with catalyst because the downstream wall temperatures are too cool to be 
effective. 
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CAPSTNE5.MCD 



cm '^0.01 m 
N '^kg m sec 



R ^8.3143 - 
K 



Gas Properties and Constants 



J^ kg m^ sec ^ 
Pa N-m ^ 



W-J-sec' 
atm = 101 300- Pa 



Pref - 1 atm 
Hm = 0.000001 m 



Rc 8.3I43 J mole **K ' 



MolWt 



'ft 

5 

ii 



n Gas 

0 CH4 

1 02 

2 N2 
3C02 
4H20 



5G7H14 

methyl-cyclohexane 



a(U) 



e(U) 



MW^j ^ 16.044 gm o^j 3.785-10* '** ni e^^ :M48.6 K 



MW, : = 31.999-gm o 3.464- 10* '^-m e, = 106.7-K 



MW ^28.0l8gm a2^3.798 10''^-m S2^71.4K 



MW3 =44.01 Ogm O3^3.94M0"***m £3 = I95.2K 

MW^ M8.0l5-gm =2.64M0'*^m =809.1 K 

Tc7 ^572.2-K Tb7 =374.1 K Pc7 ^3.47-10*Pa 
I 

Tc7 



1.2593 



-Vc7^10"*°ni 



MW5 ^98.I89-gm = 7.1661 • 10 **^ =454.3794* K 

6 CO MW - 28.01 gm - 3.69- 10'^** m E^ --91.7'K 



ru 

□ 



Define constants for transport integrals: 





= 1.16145 




^0.14874 


aDa^ 


^1.06036 


cDb, 


= 0.1561 


ajiaj 


-0.52487 




^0.7732 


aDaj 


"0.193 


aDbj 


= 0.47635 




^2.16178 




-^2.43787 


aDa^ 


M.03587 


aDbj 


= 1.52996 










aDa^ 


^1.76474 


oDb3 


= 3.89411 



Define transport integrals: 

Qm(T) -anaQ-expj - ianbQ'ln(T); ; 

nD(T) oDaQ-exp; ^oDb^-lnCT). 



2 

iM= 1 

3 



ajia.^-expi qib.^'T: 



oDa. exp oDb. T 



fifi(3) = 1.0394 



pD(3) =0.95002 



Vc7 -368 



CAPSTNES MGD - 



*««*******«**«**•*-****** ******«4 



Set Inlet Conditions: 
♦ 0.1 

0.205 



ycqc7 -f 



14 



yeqc7 =0.001464 



P ' I -atm 
yadc7 0.0021859 

, . , yadc7 

9eq -9 - 

yeqc7 

(t>eq =0.149281 



This valu e for yadc7 was determined via enthalpy 

functions^mHUHIHHHIIfoi' temperature 
rise from 1 0OOF to 1 51 1 F. 



The value for yeq was determine for a given equitvalence 
ratio of 0.1 in moist air 



Define gas composition with C7H14 fuel using moist (2%) air: 



;(l - y0) 0.205- lO.SyO x 
( 1 - yO) 0.775 
7-(y0.x) 
i (1 -y0)0.02-h7-y0-x 
' yO.(l-x) 
0 



Wy7(y0,x) 



yyeq =y7(yeqc7.1) 
pyad -y7(yadc7,l) 

IB 



1 ^ 2.5 yO-x 



yeq - 



0 

0.188634 

0.771043 

0.010213 

0.03011 

0 

0 



yad = 



0 

0.180613 

0.769103 

0.015218 

0.035066 

0 

0 



Define functions for evaluatton of mixture properties (ideal gas) for complete fuel conversion and a given temperature. 
Density and average molecular weight: 



pmix(T) 



6 

kk = 0 



_p_ 

RT 



MWmix 



kk = 0 



pmix{300 K) = 1.1 62787' kg- m ' 



MWmix =28.631 -gm 



CAPSTNE5.MGD- - 



Define drffusion coefficients for specific binary diffusion pairs: 



, ...L. .T3; 



:0.5 



Dif[(n,m.T) ^ 1.8829- 10 



i \ MW MW 

27 [ \ n n 



•10 -gm 'cm -sec -K 



■ a o \ 



\ 2 1 



.QD 



\ n m/ 



0.5 



P 

Pref 



Dif!(0,2,300-K)=0.2225-sec"* W 



Difl[(2.6.300 K)-0.207-sec ' W 



Define functions for use with various fuet-air mixtures: 
1 



QDch4(T) 



4 

V 

jgas= 1 



Difi(0,jgas,T) 



Dch4(298-K) =2.1912-10 ^ -m^-sec ' 



SJDc7hl4(T) 

01 



E 



yeq 



jgas= 1 



Difi(5,jgas.T) 



Dc7hl4(298 K) =5.8372-10 ^ 'm^'sec * 



OJdco(T) 



ry 
Q 



4 



yeq 



. ^ Dif!(6Jgas.T) 
"1 = 1 



Dco(298 K) =2.0451-10""^ -m^-sec ' 



Display gas mixture properties 



idsp 0..9 



Tin '=1 ( 1000 - 32) - ^ 273.15 K Tin = 810.9278-K 



Tdsp ^{750^ 50idsp)K Tout=|(15ll 32)- + 273.15 j-K Tout = 1094.8 167-K 



CAPSTNE5.MGD- 



4 



2.5*10 



2M0 * 



Dc7hI4{Tdsp.^^i 



Dco(Tdsp., 



no 



5'10 




1200 



idsp_ 



L M7.5'10^'K 
91 810^K 



^9.5I0^K 



m MO-K 



1.082- lO'^-m^-sec~* 


3.1097- 10' ^-m^-sec'' 


12064- 10"'*m^-sec'' 


3,4776- 10" ^-m^-sec"^ 


1.3359- lO'^m^sec"' 


3.8613-10"^-m^-sec"* 


1.4705-10''*-m^sec"* 


4.2603- 10" ^-m^-sec'' 


1.61-10 "^ m^ sec * 


4.6743- 10" ^-m^-sec"' 


1.7543 I0''''m2*sec"* 


5.1031-10"'-m^-sec'^ 


1.9033- lO'^m^sec* 


5.5462- lO'^-m^-sec'* 


2.057- 10''*-m^sec'* 


6.0035- lO'^-m^-sec"* 


2.2152- lO'^-m^-sec* 


6.4746- 10"^-ra^-sec"' 


2.378- 10"^-m^-sec"' 


6.9593- lO'^-m^-sec' 



Ml^idspL 



9.8836-10'^ m^-sec' 



1. 101 



1.2184 



1.3402 



1.4665 



1.5972 



!.7322 



1.8714 



2.0148 



2.1624 



10 



r4 



10 



r4 



10 



r4 



10" 



10' 



10 



r4 



10 



r4 



10' 



10" 



m^sec' 



m^sec" 



m^-sec~ 



m^sec" 



m^sec" 



m^'sec" 



m^sec' 



m^sec" 



m^sec" 



Define interpolation fonnulas for Dch4. Dc7h14. and Deo vs. T(K) via Jandel Sdentific's Tablecurve program 



Dch4T(T) -exp|- 10.623294 . 1.6259688-ln|l] ~ cm^ 



sec 



Dch4(1000'K) - 1.75426-sec'' W 
Dch4T(1000-K) = 1.75432 -sec"' xm^ 



Dc7hl4T(T) -exp[ IL822I89' !.6260l38.1n| -) - — ^^^^ !-^- 



sec 



Dc7hl4(l000-K) = 0.5103 Tsec ' W 
Dc7hl4T(1000-K) =0.5IO3l'sec * W 



DcoT(T) =exp| 10,831955 i 1 .6394196-lni ™ j 24.53899-- j cm^ sec ' 



Dco{IOOO-K) = 1.59722-sec ' W 
DcoT(lOOO-K) - l.597l4'sec * W 



■eAPS"fNE5.MCD— ^ 




Define viscosity functions for nriajor constituents: 



,-26 (MWoT^°^^ 10^7gm\°* 



KT) 



2.6693- 10 



.0.5 





T\ 








1"' 












,0.5 



; .2 ^ /T\ sec 



I.II6(MW^-T) 



2.6693-10 



2.6693- 10 



26/ 10^ cm \ /gm>^^ 
sec 



26f l0^-cm \ /gm^Q-^ 
sec 



0.5 



01 

w 

a 4 



(3.1l510"*cni)^-£i|i 



514-k) 



2.6693.10-".(^^j- 

sec / 



gm 
K 



0.5 



m 

3 S 



n= 1 




K300 K)- 



yeq 



1.10281-10 
2.0603' lO"^ 
1.76981-10' 
1.51851M0 
1.06784-10" 



0 

0.1886 

0.771 

0,0102 

0.0301 

0 

0 



-5 



"5 



kg'm ^'sec 



^lIl^i3l^TdspJJ ; 4-10 




Tdsp 



1200 



Tdsp 



750 K 



800K 



850 K 



900K 



950 K 



1 000-K 
i050^K 
[ilOOK 
IllSOK 
:i200-K 



^imix(Tdsp.. ) 


3.3875- 10" ^kgm' 


^•sec"' 


3.5348- lO'^kgrn" 


•sec' 


3.6787- lO^kgm' 


^-sec' 


3.8194- 10' ^-kg-m' 


*-sec'* 


3.9573- 10' ^-kgm 


•sec' 


4.0926- lO^-kgm 


*-sec"* 


4.2254-10 ^-kg m 


•sec' 




'sec'^ 


4.4844-10 ^ kg m 


**sec * 


'4.6108- 10 ^ kg-m 


'•sec'' 



idsp 



CAPSTNE5MCD 



HmixT(T) ^expf 0.6984944 r 0.640549 ln^^ j ^ li"*-^-?-^-^^^^ 



Mx(lOOO K) =4.0926-10 ^ 'kg-m * -sec"' 



HmixT(IOOO K) =4.09256-10'' 'kg-m'' -sec ' 



Define thermal conductivity functions for light molecules 

igas -0, 1 ..4 





-1.869 




0.08727 ' 




1.179- 10"* 




-(3.614-10'*) 




- .3273 




0.09966 




-(3.743 10"^) 




9,732- 10'' 




0.3919 


hi - 


0.09816 




-(5.06710'^) 


(ft- 


1.504- 10 * 




-7.215 




0.08015 




5.477- 10"* 




-(1.053-10'*) 


3 S 


7,341 . 




0.01013^ 




1.801-10"^ 




_ -(9.110 *) 



CH4 

02 
N2 

C02 

H20 



s ! 3 2 ' 



(I 



0.00001- 



w 

cm-K 



m 



1=1 



yeq.-Mi,T) 



1 



MW. 



[ / MW.\ 
8-1 1 -^ 



I 



MW., 
J/ J 



0.5 



X(igas,IO00K) 


0.16716-kg-msec 




0.07163-kg-msec' 


3.Kl 


0.06292-kgm-sec' 




0.06788 kg-m sec 


3.K-1 


0.08631 kg- m-sec' 





Xinix(1000-K)= 0.000653 



W 
cmK 



0,075 



0.065 - 




0.055 - 



0.05, 



'600 



800 



1000 



1200 



Tdsp 



7.5-10^-K 




810^-K 




8.5-10^-K 




9.10^-K 




9.510^-K 




I-IO^-K 




1.05- lO^K 




I.MO^K 




1.1510^-K 




l.2ltfK 





Xmix(Tdsp.^^)_ 



0.05323 -kg-m-sec'^'K 



0.055775 -kg-mse c'^-K 



0.058242-kg-tn-sec'^-K' 



0.06064 Ikg-m-sec^-K 



0.06298-kg-m-sec^K 



0.06 5268-kg-msec'^-K 



0.0675 12- kgmsec'^-K 



0.06972 Ikgm-sec'^-K' 
0.071904 kg-m-sec'^-K 



0.074069-kg-m-sec -^-K 



Tdsp, 



tdsp 



CAPSTNE5.MGD- 



XmixT(T) -expl 7.073383 r 0.6397005 Ini - j ^ ■ ' — j W m '-K * 



Xmix( lOOO K) = 0.0653- m ' K ' 'W 



XmixT( lOOO K) = 0.0659-m 'K ' ' W 



Define heat capactties: 



Cpcoef - 



19.25 .05213 0.00001197 0.00000001132 

28.11 -0.000368 0.00001746 -0.00000001065 

31.15 -0.01357 0.0000268 -0.00000001168 

19.8 0.07344 -0.00005602 0.00000001715 

32.24 0.001924 0.00001055 -0.000000003596 

76. 1 9 0.7867 - 0.0004204 0.000000075 1 6 



"1 



hjij =0..3 

y 

plCp(0,300 K) =35.6607-K ' J 

m 4 

' " I 

^Cpmix(T) - Y, Cp(jlc.T) yeqj^ 
^ 4 

fU J 
p3 Cpclmix^ -- XI ^P'°*'kjjy*^jkK 

n jk= 1 



jj = 0 



Cp(5,300 K) = 124.0133-K 'J 



Cpinix(998-K)= 33.5144 "J-K 



-I 



J/mol/K 



Cpcfinix = 



30.493459 
-0.009725 

2.370305- lO" 
I -1.094786-10' 



kg-m^'sec"^'! 



CpmixT(T):- Cpcftnix^-j^j^ 



Cpmix(lOOO K) =33.52413-kg'm^-scc"^"K * 
CpmixT(lOOO K) = 33.524 13" kg- 'sec"^ 'K * 



igas 0,1.. 6 



P = 1.013*10^ 'Pa 



?♦•••♦*•♦•♦♦♦♦**«♦♦*♦**»**♦*«♦♦***••***♦♦*•**♦♦♦•*♦**♦** 
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Define momentum and thermal drffusivfties: 



a(igas.T) '-^ 



X(igas.T).R.T 



Cp(igas.T)P 



a(0.300 K) =0.2344-sec * xm^ 



amix(T) 



umix(T) 



umixT(T) 



^ A.mixT(T)R-T 
CpmixT(T)P 

{^ixT(T) 
pmix(T) 

, MmixT(T) 
pmix(T) 



amix( lOOO K) = 1.6 132' sec * W 



umix( lOOO K) = 1.1732-sec ' W 



pmix( lOOO K) = 3.48836" 10 ^ 'gm cm*^ 



pmixT(T) -pinix(T) 



Dch4T( lOOOK) = 1.7543 "sec ' 'cm^ 



Dc7hl4T(1000 K) =0.5103-sec * W 



^Aefine Prandt, Schmidt, and Lewis numbers 



■in 



SI 



UJ amix(T) 



;j5cIT(T): 



omix(T) 
Dch4T(T) 



Pr(lOOOK) =0.7272 



SclT( lOOOK) =0.6688 Sc7T(T) 



_ umix(T) 
"Dc7hl4T(T) 



Sc7T(1000 K) =2.299 



, umix(T) 



^ DcoT(T) 



SccoT(lOOOK) =0.7346 



lU 
p 



SclT(Tdsp.^ 
ScTT 



0.9 



\, 0.8 



_ 2.5 
SccoT(Tdspjj5j^) 



0.7 



0.6. 



700 



800 



900 



tooo 



tlOO 



1200 



These relations have weak 
temperature dependence and 
(presumably) very weak 
dependence on composttion 



Tdsp, 



idsp 



This ends the extended section that defines gas properties for the combustion mixtures. 
*«•***•••***♦•«•**♦♦**♦*♦*♦*♦*********************♦♦***•******♦***♦♦♦**••*•**♦**♦•*♦*•*♦*********** 
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Honeycombs: Represent as square channel metal alloy wall 0.003 mil. uncoated; 0.004 mil coated 



dwalluc -0.003 in 



dwalluc =7.62-10 'm 



Cells! -lOO-in^ dwalli ^dwallc 

dcelll =(CellsI) ^ - dwalll 

dcelll =0.2438-cm ksl -0.002-cm 

Ecein ^dcelll^ Cells! eceill =0.9216 



dwallc 0.004- in 



dwallc = 1.016- 10 -m 



Cells2 -20010^ dwall2 - dwallc 

dcell2 (Cells2) dwall2 
dcell2 =0.1694*cm ks2 -0.002 cm 

£ceI12 ^ dcell2^ Celis2 ecell2 = 0.89006 



Cells4 =400in" 



dwall4 = dwallc 



dcell4 : = (Cells4)'*^^ - dwall4 dcelM =0.1 168 *cm 
dcell4 =0,I168-cm ks4- 0.001 -cm 

£ceII4 ^ dceI14^-Cells4 Ecell4 = 0.8464 



Cells6 ^600- in' 



dwall6 = dwalluc 



dceU6 = (Cells6)'°^ - dwall6 dcen6 = 0.0961 "cm 
dcell6 =0.0961 'cm ks6 ^O cm 

Ecell6 = dceI16^-Cells6 £cell6 = 0.85843 



Cells8 =800-in 



dwallS - dwalluc 



J dcell8 =(Cells8)"°^ - dwalW dcellS =0.0822 'cm 

dcell8 =0.0822 -cm ks8 -0-cm 
W£cell8 = dcell8^Cells8 Eceil8 = 0.83749 



01* 



r****** ••♦♦♦♦♦***♦**♦♦*»**♦♦♦♦*♦**♦*♦♦♦** 



g1 Define recuperator and monolith configurations: 

dreco =17 in 



dreci =6- in 
035d M.375-in 
dhco - dreco - fid 
dhci = dreci ^ od 



dreco =43.1 8 'cm Lrec =7in 

dreci = 15.24 "cm Lrec = 17.78 'cm 
Add space for mounting and insulation? 

dhco = 39.6875 'cm Lhc : ^ 4 in 

dhci = 18.7325 'cm Lhc = 10.16 'cm 



She 



dhco^ - dhci^ 



She = 0.0961 "m 



Srecup =0.0958 m^ 



Define initial combustion conditions 



jcase -0,1.. 4 



5M= 0.9964 
She 

♦**♦♦♦♦♦««*********•♦♦••*«* 



0.2 1350 340 I 

\ 0.26 1350 390 ! 

case ! 0.314 1350 440 1 

0.38 1350 490 ; 

= 0.72 1190 675 I 



Mass flow (lb/sec) 
Recup T inlet (F) 
Recup T outlet (F) 



TKF(T) - 273.151; r 32 
\K ; 5 



TFK(.) (•^f■5•^ 273..5.K 
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Prec ■ I atm 



Phc - 1 atm 



mf , 



case. „ 0.454 kg sec ' 

lease, 0 ^ 



Trecin. . , -TFKfcase. . i 

jcase-M V jcase.l/ 



Trecout. , , ^TPKicase. .) 

J case +■ I \ jcasc.2/ 



0 I 

0.0908 I 

0.118 ! 

I' kg* sec 

0.1426 
0.1725 
0.3269 



Urec. 



mf 



Trecin = 



0 

1005.37 
1005.37 
100537 
1005.37 
. 916.48 



jcase-H t 



pmixi Trecin. ^ Srecup 



'^nly half of the recuperator cross section 
'^is available for exhuast gas flow 



Urec 



Trecout = 



444.26 
472.04 
499.82 
527.59 
630.37 



0 

5.4633 

7.1023 

8.5774 

10.3803 

17.929 



'm'sec 



Grec 



mf 



jcase + 1 



jcase -H 1 

Srecup^ 

2 / 



9^ 



mf . , 
jcase -t- I 



jcasc -t \ 



pmixfTrecin. ,VShc-6ceH2 



Uhc4 



jcasc -|- 1 



■■ Uhc2 



jcase t I 



6cell2 
ecell4 



Uhc6. 



jcase + 1 



= Uhc2 



jcase -J- 1 



6cell2 
ccell6 



!lJ 

&hc2 = 



0 

3.058 

3.9753 

4.801 

5.8101 

10.0353 



•m'sec 



Uhc4 = 



0 

3.2157 
4.1804 
5.0487 
6.1098 
10.553 



' m'sec 



Tavg - 



Trecin ^- Trecout 



I 0 



724.8167 



Tavg = 



738.7056 
752.5944 
I 766.4833 
! 773.4278 



Define average temperature etc/ for 
estimated AP and AT in the recuperator 



pmixreciTavgj j = 0.48 1 3* kg* m 
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******************** 



Evaluate Reynolds numbers for use In calculating the momentum, heat, and mass transfer functions: 



Rercc6(G,T) 



MmixT(T) 



Rerec6(GreCj,Tavg j =54.9906 



Rerec8(G.T) 



GdcellS 
MmixT(T) 



Rerec8(GreCj.Tavg j =47.0389 



For laminar flow with entrance effect (adapted from Kays & Crawford, "Convective Heat and Mass Transfer^: 



clTicrec8(G,T,z) 



14.227 



Rerec8(G,T) 



I -h0.046263-|Rerec8(G,T)-^^j 



0.45363 



pfricrec6(G,T.z) - 



14.227 



Rercc6(G.T) 



I 0.046263 |Rerec6(G,T)~'^i 



Average coefRctent of friction 
0 45363 a smooth square channel 



; Calculate monolrth pressure drop in recuperator with entrance effect: 



(Urec. . ,) 
\ lease + I / 



ii&Prec8. , , -4 cfricrec8/Grec. , i,Tavg. ^,,Lrec) 
^ jcasc + 1 V jcasc -i- 1 * ^jcasc -hi' / 



pmixrec(Trecin.^^^j) ^ pmixrec(Trecout.^^^^,)* 



Treci 



Tree 



m 
m 



(Urec. 

\ jcasc i- 1 / 

2 



/Treci 

pmixrec(Trecin.^^^ ^ ,) pmixrec(Trecout.^^^ ^ i)' [l^ 



APrecS 
Free 



0 i 
0.0097 : 
0,013 I 
0.0162 [ 
0.0203 I 
i 0.0393 ; 



APrec6 



0 

0.0071 
_^ 0.0095 
Prec I 0.0119 
i 0.0149 
I 0.0289 



APobs P 



! 0 

jO.Oil 
i 0.014 

10.018 

! 

i 0.022 
; 0.042 



Chose 800 cells per in2 
with 10% length factor 
in Lrec for chevrons 
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Averaged Nu for 800 cpi square channels (fully developed 
laminar flow) taken from Rosner: 



7.6 (Lrec Ll) 



Rerec8(Grec. . ..Tavg. .,)Pr(Tavg. . .VdcellS ! 

V jcasc -t- 1 * *'jcasc -k- 1 / V ^jcase 4-1/ j 



NuhavgrecS 



^Jirecavg, 

E . i ' 

u 



0 

2.976 
2,976 
2.976 
2.976 
2.9762 



Entrance effect for heat transfer within the recuperator is only 
significant for the first half inch - average Nusselt No. 2.976 for all cases 



Xinix(Tavg.^^ ^ j ) NuhavgrecS^ 



jcase -t- 1 



*jcase +• I 



dcellS 



hrecavg = 



0 

188.0022 
190.6356 
193.2448 
195.8307 
197.1273 



'kg'sec * 



WMi2/K 



WATskin 

m 

fW 

m 



Grec, •dceii8(Cpniix(Trecout. ^ , )Trecout 

jcase -f \ \ \ jcase -t- 1 / j 



jcasc -h 1 



CpmixfTrecin , , )Trecin. .i 

*^ V jcase y \l jcase j- X I 



jcasc 1 



□ 



ATskin = 



0 

-7.5657 

-9.274 

-10.5339 

-11.9549 

-13.635 



MWmix hrecavg^^^^ ^ j*4 Lrec- 1.1 



0 

-27.2364 
_ -33.3864 
5K -37.9221 
-43.0378 
-49.0859 



ATskin-2- 



Approximate average AT (F) 
across wall In recuperator 



AT seenrts low - more conservative calculation uses dcell for 600 cells^n2 - in this case AT maximum 
Increases from 50F to about 70F 



j 0 I 

I -37.2231 1 

^ 9 /dcelI6;^ | 45.628 | 
ATskin-2 1— ! =! 



5 K \dcelI8/ 



51.8268 
"58.8182 
67.0839 



/ dee 1 16 \ 
''dcellS 



1.3667 
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************* 



Use 600 cells/in2 in recuperator dimension for conservative calculation of CH4, C7H14, and CO conversions 



Find mass transfer rates via Nusselt formulas: 



Numavgrec6cl. , , -2.976 



^ i 7.6 Lrec l.l \ 



Numavgrec6cl = 



Numavgrec6c7.^^^ " 2.976 



7.6-Lrec-I.l 



Rerec6/Grec. ^,,Tavg. ^ , )-Sc7T(Tavg. .j-dcen6 
\ \ jcasc -H 1 * "jcase \ ^jcase -h U 



Numavgrec6c7 = 



0 

2.976 

2.976 

2.9761 

2.9761 

i 2.9765 
' 0 

2.9765 

2.977 

2.9775 

2.9784 

2.9891 



y 

ni4umavgrec6co.^^ ^ i - 2.976- 

w 



1 ^ 



7.6Lrec-l.l 



Numavgrec6co = 



^pnce again the entrance effects within the recuperator are negligible 



0 

2.976 
2.976 
2.9761 
2.9761 
[ 2.9766 J 



\ y 

p fThe average OntegraO Nusselt numbers at 1-atm show very little contribution of the thermal entry length, i.e. the flow 
^velocities and heat and mass transfer coefficients have values close to those in long channels. 

Now detemnlne the maximum possible CH4, C7H14, and CO conversions in a 600 cells per in2 recuperator 



ych40 ^ 0.005 yc7h 1 40 ' 0.0005 ycoO = 0.005 

- 4Lrec- 1.1 •Numavgrec6cl 



ych4rcc6 



jcase -h 1 



^, ych40-exp r 

\Rerec6(Grec.„^^,,Tavg.^^^,).SclT(Tavg.^^^,).dcell6 



Start with maximum 
concentrations of added 
pollutants 



4LreC' 1 .1 •Numavgrec6cl 



jcase 1 



Rerec6iGrec. . .,Tavg. ^J-SclTlTavg , i dcell6 

\ jcase I r ='jcase + 1 / \ ^jcase 4 I / 



65^4083 
50.9771^ 
4275 52 i 
35!7764; 
19.0017- 



ych4rec6 
ych40 



i 0 

!o 

i 

0 
0 
0 

5.5932- 10 
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yc7hl4rec6. 



-yc7hl40expj 



- 4-Lrec-l.l-Numavgrec6c7. 



^""^^^ ' Rerec6(Grec.„. . ,,Tavg. ^ , jScTT Tavg. .jdcell6j 



J case H 



jcase ' 1 

'jcasc 1 



yc7hl4rec6 = j 



3.5287- 10 
2.1834- 10 
2.2859-10' 
1.6805-10" 
2.0597-10" 



12 



-10 



yc7hl4rec6 
yc7hl40 



0 

7.057-10 
4.367' 10~ 
4.572-10" 
3.361-10 
4.1I9'10' 



-9 



-5 



4-Lrec- 1 .1 •Nuniavgrec6c7 



jcase 1- 1 



Rerec6(Grec.^3^ . r Tavg.^^^ ^ , ) •Sc7T(Tavg.^^^ ^ ,) .dcell6 



18.7692 



14.6441 



12.2956 



10.3007 



5.492 



S ycorec6. 



jcase 1- 1 



'- ycoO-exp 



4Lrec- 1.1 •Numavgrec6co 



jcase -h 1 



Rerec6[Grec. ^,,Tavg. ^ ,)-SccoT(Tavg. ^,)-dcelI6 
\ jcase -h 1 * °jcase -i- 1 / \ ^jcase +■ I / 



01 

m 

b^corec6 = 



0 
0 

0 
0 
0 

1.457' lO' 



to 



ycorec6 
ycoO 



0 
0 
0 
0 

6.4676-10' 
2.914*10"^ 



15 



* 4-Lrec- 1 .1 •Numavgrec6co 



jcase + 1 



Rerec6(Grec. ,,Tavg. ^ . )-SccoT(Tavg. . , Vdceil6 

\ jcase -t- 1 ' ^'jcasc -I- 1 / \ ^jcase -f- 1 / 



59.7854 



46.5803 



-39.0561 



32.672 



17.3512 



o 

iJA/ith SOOOppm C7H14 in turbine exhaust, exit concentration reaches 20 ppm at the highest power level 
' for conservative 600 cells/ln2 calculation. For 800 cells/in2, the exit concentration in this case would reach 
about 2 ppm. With only SOOppm C7H14 In the inlet, should get <2 ppm in the outlet of the recuperator 
assuming favorable catalytic rates at the wall. Also this calculation used averaged temperature and flow 
conditions - will compare below a 1 -dimensional (axiaO calculation tiiat determines the effect of decreasing 
temperature on the kinetics and mass transport. 



Maximum possible (transport-limited) CO and CH4 conversions in a 7 in. 
600-800 cells/in2 recuperator are >99.99%. while cyclohexane conversion is 
>99.6%. The final CO and C7H14 levels are < 1 ppb and .2 ppm for post 
turbine concenti'atipns of 1000 and 50 ppm. respectively. 
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An alternative method of calculating perfornnance uses local Nusselt nunnbers for heat and mass transfer to determine the 
wall temperature and combustion rates at the wall. First we define the local Nu function then check conversion assuming 
fast kinetics and compare with the previous averaged (integrated) Nu formula. 



Local Nu(x) taken from Kays & Crawford for hydro & thermal laminar flow entrance 
with circular channel (modified for square channel by a ^ctor 2.96/3.66: 



ncase ~ 5 



Numono(2) -(I t exp(- 1.09689 - 0.38140Mn(z) - 8.287774 z**^))-2.976 Numono( 0.1 152163) =3.112 



zzchar6c7 



Rerec6(Grec. ^..Trecin. ,) Sc7T(Trecin. . , ) dcell6 i™ 

. \ J case -t- 1 jcase -t- 1 / \ jcase j- 1/ L>fCK^ 



jcase f I 



zzchar6c7 



= 0.9444 



Case specific 



Rerec6(Grec. ^,,Trecin. ^, )'SclT Trecin. ^ , •dceII6 i 

zzchar6cl - ^ ^-^^ ^^^^ — 3.2574 

jcase ^ I 2 zzchar6cl 

ncase 

finz =400 
j^^iz =0..nz 

^'*=' Lrcc 

m zrec6. - iz ^ dwall6 

« nz 



Lrec = 1 7.78 'cm dcell6 ^ 0.0961 -cm 



zrec6. 



zrec6. 



yj 
W 



zz6c7. =- 



*^ zzchar6c7 



zz6cl. " 



^ zzchar6cl 



zzchar6c7 - 



Case specific 



0 ! 

4.896 

6.3647 

7.6867 

9.3023 

18.8265 



□ 
i fa? 



NuiTiono^zz6c7-^ 
Numono (zz6c I ■_) 




no 



zzchar6cl 



0 

1.4244 
1.8518 
2.2364 
2.7064 
5.4583 



cm 



Case specific ncase = 5 
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Case specific ncase = 5 

Dc7hI4T(T)Numonoi - 

* zzchar6c7 

kgpbl6c7(T,2) ■ """^ kgpbI6c7(Trecin .Lrec) = 13.61 59-sec ' 'cm 

dcell6 "case / 

Dch4T(T)Nufnonoi " ~~ t 

1 2zchar6c 1 ; 

kgpbI6cI(T.2) ■ — "^^^ kgpbI6cl{Trecm ,Lrec) =46.9581 'sec * -cm 

dcell6 ^ ' 



Trecin - Trecout 

Trecz - Trecin !^^.iz Trecz„, = 630.3722-K Trecout = 630.3722'K 

iz ncase ^ nz ncase 



Determine maximum conversion of fuels with very fast catatytic combustion rates (zero fuel concentration at the wall) 



^zc7hl4 -yc7hI40 



yzch4^ = ych40 



^yzc7hI4^^j ^ jl - kgpbi6c7(Trec2.^.zrec6.J-— 1.1 ^ 



: t I 



ecell6 Trecin^^^ 



nz dcell6 Urec Trecz. / 

ncase iz / 



•yzc7hl4. 



ncase = 5 



Case specific 



y yzch4.^^j :^ j 1 - kgpbI6cl(Trecz.^,zrec6j-~ l.l- 



ECeIl6 '^^^^%case 



nz dcell6 Urec Trecz. 

ncase 12 



•yzch4. 



Case specific 



H'zc7hl4^ = 3.385- 10 " yzcThM^ =0.0005 yzc7hl4 
^ "'=6.7710~' 



nJ y2ch4^ = 2.1646- lo '" yzch4Q = 0.005 



yzc7hl4. 



yzc7hI4Q 

yzcM.^ 
yzch4„ 



0.01 



no 



no 



yzc7hl4 




Case specific ncase = 5 



yzch4^ 
yzch4„ 



= 4.3292-10 



Conversion of C7H14 under case 5 
(full power) is 99.3% [vs. 99.6% via 
averaged integrated formulas] with 
about 1 inch entry length. Conversion 
of CH4 (and also CO) is very high 
(>99.99%) within 3.5 inches. 



This result shows that the averaged 
integrated values for CH4, C7H14. and 
CO conversions are very good 
approximations and can be used for 
drfferent values of fuel and CO 
concentrations. 
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Add finite rate kinetics for PdO/AI203 methane combustion catalyst 
Begin by defining a few constants and catalyst functions 

R-82.076 K~* -atm cm^ atm*cm3/morK 

Rem ^83l43000 gm cm^'sec ^ K * TOL ^O.OOOI 



Define Catalyst Functions: 



nsph(*) -^.(— L--- ij 



♦ \tanh«») 
tanh(^) 



3 ♦ 



i Catalyst Physical Properties for Pd/AI203 



y'Scat -85 m^ gm'' pcatpart =2.5 gm cm'"' 

^Mpore :-0.I5-cm^ gm"' dpcat : = 0.0002-cm 

gxat= 0.0020 cm 

nj 

!w , Vupore 

?^pore: = 6 — efipore =0.375 x^ipore -3.0 

— Scat 

cftipore== 1.0588- 10 ^ 'cm 



Def!CH4(T) ^ - 



gppore 
Y^pore 



I 



Dch4T(T) 



dupore- 



•RT 



#pore' 



•R-Tin 



= 0.l095-sec * W 



DeffCH4( Tin) =4.43271-10 ^ sec* W 



Dch4T(Tin) = 1.2343 -sec ' 'cm^ 
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Def!C7H14(T) 



ejipore 
ippore 



j Dc7hl4T(T) 



dMpore- 



•RT 



j^.RTin 

dupore- \~ =0.0829-sec * 'cm^ 



Def!C7H I4( Tin) =3.20528-10 ^ 'sec ' W 



Dc7hl4{Tin) =0.356"sec ' W 



DmacCH4(T) :-^ Dch4T(T) 
DmacC7HI4(T) •Dc7hl4T(T) 



J- kOT = LO gm * sec ' 



DmacCH4(800-K) =0.193'sec W 



DmacC7H14(800-K)-0.0556-sec ' 'cm 



Test value for kOT 



ui Microscopic Thiele Modulus for CH4 

01 Mk.T) ^ ( pcatpart 



/R'T'\ k 

\Pref/ DeffCH4(T) 



dpcat 
6 



<t»KkOT,Tin) =0.2042 



Q 

Q Macroscopic Thiele Modulus for CH4 

^ 4>mac(k,T) ^ 

□ 



pcatpart.l^^j.k 



T]sph(»Klc>T)) 



DmacCH4(T) 
Microscopic Thiele Modulus for C7H14 



Lcat 



♦n7(k.T) vl[pcatparf(Mj 



Def!C7HI4(T) 



dpcat 
6 



<I>mac(kOT, Tin) = 1.8333 



♦n7(kOT,Tin) =0.2402 



Macroscopic Thiele Modulus for C7H14 

«l>mac7(k,T) 



peatparf(M:.,|. nsph(»ti7(lc.T)) ^^^^ 

'Pref : DmacC7H14(T) 



I pcatpart- j- 



a>mac7(k0T,Tin) =3.4117 



RatekCH4(k,T) ^nsph(Wk,T)) nfpl(<l>mac(k,T))-i pcatpart | --! k 

; \Pref/ ] 



RatekC7H14(k,T) ^ nsph(«^7(k,T))-nfpi(Omac7(k,T)) 



R-T 

pcatpart- j -k 

\Pref; 



RatekCH4( kOT , T in ) = 85999.49 1 4 * sec 



RatekC7H14(kOT,Tin) =48479.296'sec 
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Defines transport restricted rate constant function given catalyst effective first order rate constant, 
kOT, methane or cycfohexane concentration, x, and surface temperature, T 

PdO Catalyst Rate Constants 

ErpdO =I2800 K Erpdl =4900-K Erpd3 M5000 K 

ReglonO - low T PdO to 625K Regioni - Intermediate T Region2 - Flat 909K Region3 - 
ignores inhibition by H20 and C02 625K to 909K to 1 048K decomposition Pd metal 



ApdO :-2.4 I0^.200 gm -^sec * Apdl - ApdO-expf-^^2^ ^- ^^^^ j 

\ 625-K 625-K/ 



Apd2: = Apdl.exp(^^2^) 



ARd3 = Apd2-exp 



Erpd3 I 
1048.1. K 



Interpolation formula for the four PdO+Pd rate constants 



^pd(T) 



: - r 



Apd3^.exp(-2:^j^ 



(Apd2-^^ Apdr='.exp(2:M2dl j ^ Ap^o.2.,^p|2:^jj 



0.5 



T :=500Ki-nx-25-K 
nx 



Onx =0,1..41 

w 

CATALYTIC RATE CONSTANT vs TEMPERATURE 
Q (mol/atmCH4/s/gcat) 



HI 




lOOOK 
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keffrec6cl(T,z) - - - - 

iRatekCH4(kpd(T),T) Lcat/' ^ kgpbl6cl(T,z) 



kgpbI6cI(Tavg^^ 



(Tavg^„„,— ! =35.3023-sec ' 'cm 



2 / 



keffirec6cl^Tavg^,^j =0.7922 
kef&cc6c7(T,z) 



sec ^ 'cm 



RaJekC7H14(10 kpd(T),T) Lcati ^ kgpbI6c7(T.z) 



RatekCH4(kpd(Tavg^^^j,Tavg^^j-Lcat =0.8l04-sec 'cm 



kgpbl6c7[Tavg^^,^y^ 1 = 10.1756-sec' ' *cm 



kef&ec6c7iTavg^^,^j =4.2778' 



sec "cm 



10 



J3 



y I 

I ! 3 



U3 



kef&cc6cl (Trccz^.zrec6j^ 

RatekCH4 (kpd (Trccz-^ . Trecz ■ j • Lcat q j 

kgpbl6c I (Trecz^. zrec6.^ 



0.01 



no 



Case specific ncase=5 



RatekC7H14(l0-kpd(Tavg^^^).Tavg^^) Lcat-7.3807'sec ' 'cm 




0.05 



0.1 
zrec6. 



0.15 



0.2 



fceffrec6c7(Trec2^,zrec6.J 
RaiekC7HI4^10 lcpd(TrccZj^ .Trcc2|^) Lcat q , 
kgpbI6c7(TrccZj^. zrcc6. J 



0.01 
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Determine conversion of fuels with transport-limited finite catalytic combustion rates: 



21 



•***••******•****«••*««*** 



Trecin - Trecout 

.... ncase ncasc . 

Trecz. -Trecin iz 

iz ncase ^2 



Trecz =630.3 722' K 

nz 



Case specific ncase = 5 



Trecout = 630.3722- K 

ncase 



yzkc7hl4. =yc7hl40 



yzkch4Q ■ - ych40 



yzkc7hl4 



kefBrec6c7 ( Trecz. ,zrec6. — 

^ »^ nz dcell6 Urec 



Trecin. 



yzkch4. , , = ! I - keffrec6cl{Trecz. ,zrec6. ) -^^ l.l — ~ 
^ iztl ^ \ IZ' nj „2 dee 



Trecz. / 
ecell6 ^^^^^s\ 



i •yzkc7hl4. 



n 



^Oyzkc7hl4 =7.036-10 



yj yzkch4^= 0.0034 



yzkc7hl4Q =0.0005 
yzkch4Q= 0.005 



dcell6 Urec 



yzkch4 

Trecz. ' 



yzkc7hl4^ 
yzkc7hl4- 



= 1.407-10 



yzkch4^ 
yzkch4 



= 0.6838 



Case specific 

ncase = 5 

Case specific 



01 



yzkcThU- 



fy yzkc7hl4Q O.I 



yzc7h!4Q 



0.01 



no 



r3 



0.05 



O.l 
zrcc6. 



0.15 



0.2 



Case specific ncase = 5 

Conversion of C7H14 under case 5 
(fiiil power) is 99.4% (Vs. 99.3% for 
the maximum BLtransport-^limited 
conversion] with gradually decreasing 
(rate limited) conversion. Conversion 
of CH4 is very poor (35%). CO 
conversion is prot>abty very high given 
higher catalytic and transport rates 
than with C7H14. 

This result shows that the specific 
rates of catalytic oxidation of all added 
components needs to be measured as 
applied to a recuperator. The use of a 
recuperator for emission control would 
appear to be greatly enhanced with 
higher exhaust (recuperator) 
temperatures. 



In the following section, the conversion of C7H14. CH4 and CO Is considered for various honeycomb monliths 
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The 



jcase t 1 



■ Trecin. 



jcase ! 1 



pmixhc. 



jcase t- 1 



■ pmixreci Trecin. 



jcase t- 1/ 



Ghc2 



'jcase i- 1 



mf 

• jcase -h 1 

ShcEcell2 



Ghc4. 



jcase 1 



mf 

Shcecell4 



Ghc6 



jcase +• 1 



mf 

: jcase + I 

Shcecell6 



Uhc2 = 



0 

3.058 
3.9753 
4.801 
j 5,8101 
[ 10.0353 



m*sec 



Uhc4 



0 

3.2157 
4.1804 
5.0487 
6.1098 
10.553 



•m'sec 



Uhc6- 



0 

3.1706 
4.1218 
4.9779 
6.0242 
10.4051 



m'sec 



•£i Evaluate Reynolds numbers for use in calculating the momentum, heat, and mass transfer functions: 

w 



01 Rehc2. 



01 

2 



jcase 1 



Ghc2. ..•dcell2 
. jcase -f- 1 

»""«T(Thc.^„^,) 



Rehc4 



Ghc4. /dcell4 

. jcase 1- 1 

j""'^* 'limixT(Thc. 

*^ \ jcase + 1/ 



m 


0 




0 




0 


lU 


43.776 




31.7427 




25.7355 


St i 

'M Rehc2 = 


56.9088 
68.7283 


Rehc4 = 


41.2655 
49,836 


Rehc6 = 


33.4562 
40.4048 




83.1744 
. 167.4548 . 




60.3111 
. 121.4242 




48.8975 
. 98.4453 



Rehc6 



jcase -t- 1 



Ghc6. , .•dcell6 

. jcase 1 

MmixT(Thc.^.^^,) 



cfirichc2(j,z) - 



14.227 
Rehc2. 



1 + 0.046263- j Rehc2.- j 



0.45363 



cfrichc4(j.z) - 



14.227 
Rehc4. 



J 



1 ^- 0.046263- 1 Rehc4 



dceII4\T-^"" 



Average coefficients of friction 
for smooth square channels 



cfrichc6(j.z) 



14.227 



Rehc6. 
J 



1 ^ 0.046263- iRehc6. 



dcell6 \ 



\ 10.45363 



J 2 / 
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^^^^^2. ^ , - 4 cfrichc2(jcase f I , Lhc) ^^-^^.pmixhc. 



APhc4. 



jcase 



^4-cfrichc4(jcase-h l.Lhc)-^ ^^^^^^^-pmixhc.^^ 



Lhc 
dcelM 



.Uhc6. 

APhc6 , : = 4 cfrichc6(jcase+- I.Lhc) ~ ^"^-^'^^ ■pmixhc 



Lhc 



APhc2 
Phc 



0 

0.0013 

0.0017 

0,002 

0.0024 

0.0041 



APhc4 
Phc 



0 

0.0028 
0.0036 
0.0044 
0,0053 
0.0088 



APhc6 
Phc 



0 

0.004 

0.0053 

0.0064 

0.0077 

0.0127 



APhcmax P - 



0 

0.003 
0.004 
0.005 
0.006 
0.011 



bJ 



Chose 400 cells per in2 
honeycombs via estimated AP 



• * ********* 



J Use 400 cell8An2 in honeycomb monolith for calculation of CH4, C7H14, and CO conversions 

m 

fU Rnd mass transfer rates via Nusselt formulas: 

□ 



Nuniavghc4cl.^^^j -2.976- 



I i- 



7.6-Lhc 



Rehc4' ^ -SclTfThc. ,)-dceH4 

\ lease f 1 \ jcase + I / t 



Numavghc4cl = 



0 

2.976 

2.976! 

2.9761 

2.9762 

2.9774 



[ 

Numavghc4co. , , ~ 2.976- j I 4- j - — r-^ :'r ':~ 

' I iRehc4 -SccoTfThc. , . dcelM ; 

I. \ lease i- 1 \ jcasc f I / 



Numavghc4co = i 



I 2.976 
' 2.9761 



; 2.9762 
2.9763 
2.9778 



•«••«••********«*«*•«***•#«******«***«••«'**•*••**•**•***••*•**-*«****•«**•#**•*«****«*«******«*•**•* 
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Numavghc4c7.^^ ^ j "2.976 



7.6-Lhc 



\ 3 



\ Rehc4. ^ .Sc7T(Thc. ^ , rdcelM 1 

\ jcase +• 1 \ jcasc i I / / 



Once again the entrance effects within monoliths are generally small 



Numavghc4c7 = 



0 

2.977 
2.9781 
2.9795 
2.9817 
I 3.0123 



Now determine the maximum possible CH4, C7H14, and CO conversions in a 400 cells per in2 monolith 



ych40 -0.005 



yc7hl40^ 0.0005 



ycoO =0.005 



1 - 4 • Lhc • Numavghc4c 1 . 

ych4hc4. . = ych40 exp ; ^-^^^ 

jcaseti *^\Rehc4. ^ -SclT THc. ^,)dcell4 

^ \ jcase -I- 1 \ lease -H I / 



4-Lhc-Numavghc4c 1 



In 



jcase -h 1 



'Rehc4. -SclTiThc 
i ' jcase ^ I 

fay 



1 l hc. , , )-dcell4 

V jcase -H I / 



ff) -48.7678 



37.5142 



-J 31.0633 
25.6688 



L;!- 12.7326 



IS 
RJ 

Pyycohc4.^^^j -ycoO-exp 



- 4-Lhc-Numavghc4co.^^ ^ ^ \ 
Rehc4. ^ , • SccoT (The. ^ , ) • dceII4 

jcase -h 1 \ jcase f I / 



ycohc4 ^ 
ycoO 



1.4726-10 
5.2301" 10 
7.0997- to' 
9.1398- 10~ 



15 



13 



ych4hc4 _ 
ych40 



0 
0 
0 

3.2315-10 
7.115110 
2,9532-10 



14 



'12 



-6 



4 • Lhc N um avghc4co. 



Rehc4. -SccoT (The. , .VdcelW 
jcase + 1 V jcase + 1 / 



44.3966 



34.1518 



-28.2792 



23.3684 



11.6029 



CAPSTNE5.MCD 



yc7hl4hc4 



jcase ' I 



yc7hI40exp 



4 • Lhc- N umavghc4c7. 



i 

yc7hl4hc4 ^1 
yc7hl40 i 



6.8527*10 
1.8058- 10~ 
0.0001 
0.0006 
0.0239 



jcase ^ I \ jcase > 1 



jcasc j- J 

)dcell4 



4Lhc*Numavghc4c7. 



jcase ^ I 



14.1935 



10.9219 



9.0478 



7.482 



3.7348 



Only for C7H14 and only In case 5 (maximum power) does the BL transport-only 
calculated conversion (97.6%) fall below 99.9% with a 400 cells/in2 monolith. With 
50 ppm C7H14 in the exhaust in this case, the exit concentration would be 1 2 ppm. 



For the final analysts, consider the finite catalytic conversion of C7I-114 and CH4 in a monolith: mease - 5 



^ zzcharhc4c7. 

jcascfl 

w 

III 

^*zzcharhc4cl. 

ffi jcase 1 1 



nz=400 



■ ^^^^^jcasc^rSc7T(Thc.^,,)dcell4 
2 

Rehc4. ^ .SclTfTrecin. ^,) dcell4 

jcase -j-l \ jcase + 1 1 

2 

Lhc = 10.16 -cm dcell4 -0.1 168 'cm 
2hc4. 



Lhc 



zzcharhc4c7 



Lhc 



zzcharhc4cl 



'0.6199 



= 2.1382 



y 
o 

0- Lhc 

n I zhc4 . =12 ^ d wal 14 zzhc4c7 . = 

^^-^ nz 2zcharhc4c7 

fi I mease 



zhc4. 



zzhc4c 1 . 



zzcharhc4cl 



NunH>no^zzhc4c7^ 
Numono{zzhc4c) 




zzcharhc4c7 



0 • 

t 

4.2621 1 
5.5407 : 
6.6914 I 
8.0979 
I 16.389 ! 



zhc4._ 
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Dc7hJ4T(T)*NumonOj' 



kgpbIhc4c7(T.z) 



\ Z2charhc4c7 / 

\ mcasc/ 



dcell4 



Case specffic mease = 5 



kgpblhc4c7(Thc^^^,Lhc) = ! I.2014-sec * 'cm 



Dch4T(T)Numono( 



kgpblhc4cl(T,z) 



I zzcharhc4c 1 

\ n 



dcell4 



kgpbIhc4cI(Thc^^.Lhcj =38.6127'sec 'cm 



keflnic4cl(T,z) 



( 



I ! 



\RatekCH4(kpd(T),T)'Lcat/ kgpblhc4cl(T.z) 



1 



kgpblhc4cl(Thc — )=38.615-sec ' 'cm 



kefflic4cl(Thc , — ) = l.8533 sec * cm 



RatekCH4(kpd(Thc ).Thc i-Lcat = l.9467*sec ' -cm 

V \ mease/' mease/ 



W kefBic4c7(T.z) =- 



ui 



I 



RatekC7H14( 10kpd(T),T)Lcat/ kgpblhc4c7(T,z) 



I 



kgpblhc4c7(Thc^^.^|?) = 



= 11. 2528 -sec 



»_kefflic4c7(Thc^^, — ) =6.7062 sec ' -cm RatekC7H14(l0kpd(Thc^^).Thc^) Lcat = 16.5977-sec ' -cm 

5 \ 2 / 



m 

rF'i 



kcl!bc4c7|'nw^jj^,zhc4.j 

RaickC7HI4i 10.kpd(Thc^j^jg) .Thc^^).Lc8t q , 
kgpbI6c7Cnic^^.zhc4J 



0.01 



r 



0.02 0.04 0.06 0.08 0.1 0.12 
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Determine conversion of fuels with transport-limited finite catalytic combustion rates: 



Case specific mease = 5 

y2khc4c7hI4Q ^yc7hl40 y2khc4ch4Q ^ych40 

yzkhc4c7hl4 =(l - kef!hc4c7(Thc _,2hc4- j.^!?—^ ???!!l_).yzlehc4c7hl4. 

\ nic«c izj „2 dcell4 Uhc4 
« mease/ 



Thc = 



i 0 
1005.3722 
1005.3722 
1005.3722 
1005.3722 
916.4833 



yzkhc4ch4. ^ , := ( I - kef!hc4cl (The .zhe4. y^^-^ |.y2khe4ch4. 



Case specific 

mease = 5 



yzkhe4c7hl4^ = 7.281- 10 



yzkhc4c7h 14^^=0.0005 



yzkhc4ch4^ = 0.003 



yzkhc4ch4 = 0.005 



yzkhc4c7hl4 



yzkhc4c7hI4 



"^-1.456*10 ' 



yzkhe4ch4^ 
yzkhc4ch4. 



= 0.5959 



01 




p3 yzkhc4c7hI4^ 
fU yzkhc4c7hl4(, 

y2kbc4cb4. 



0.1 



yzkbc4cb4Q 



0.01 



no ^ 



0.02 



0.04 



0.06 
zhc4. 



0.08 



0.1 



0.12 



Case specific mease = 5 

Conversion of C7m4 under case 5 
(full power) is 99.4% [vs. 99.6% via 
averaged integrated formulas] with 
about 1-inch entry length. Conversion 
of CH4 (and also CO) is very high 
(>99.99%) within about 3.5-Inch. 

This result shows that the averaged 
integrated values for CH4, C7IH14, and 
CO conversions are very good 
approximations and can be used for 
different values of fuel and CO 
concentrations. 



CAPSTNE6.MCD ^jj^^^^^HH 

Defines transport restricted rate constant function given catalyst effective first order rate constant. 
kOT. methane or cyclohexane concentration, x. and surface tennperature. T 

PdO Catalyst Rate Constants 

ErpdO - 12800 K Erpdl -4900 K Erpd3 - ISOOO K 
RegionO - low T PdO to 62SK Regloni - Intermediate T Region2 - Flat 909K 



ignores inhibition by H20 and C02 625K to 909K 



Regions - 

to 1 048K decomposition Pd metal 



ApdO ^2.4-10^-200gin"*-sec'* Apdl = ApdO-exp(-^^^ f Erpdl^j ^^^^ Apdl-exp 

\ 625'K 625'K/ 



Interpolation formula for the four PdO^-Pd rate constants 

kpd(T) - 



Erpd! 
909.1 K 

Apd3 = Apd2'exp 



/ Erpd3 
1 1048.1 K 



Apd3^-exp[-^). 



I 



0.5 



Apd2-^. Apdr^.exp(2:l^) . ApdO-^.exp(?:^)j 



T / ' ' \ T 

Add expression for R(.5%}Pd(3.0%) catalyst with C3H8 here, use x2 for estimated rate constant with C7H14: 
f Rate constant via crushed powder results for C3H8 



'ff, kptpdC3(T) :=6.9.10* exp(- ■^j g^ ' sec ' 



10kpd(900K) =49.2392-kg '-sec ' 



2kptpdC3(900K)=886.4358-kg '-sec ' 



5iix^0,L.41 



T :=500K + nx-25K 
nx 



m 
ru 

m 



CATALYTIC RATE CONSTANT vs TEMPERATURE 
(moVatmCHx/8/gcat) 



.O.kpd(Tj 
2.kp.pdC3(Tj 




no 



CAPSTNE6.MCD 



keffrec6cl(T.z) " — 



UateIcCH4(kpd{T).T) Lcat) ^ kgpbl6cl(T,z) 



kgpbl6cl|Tavg^^^^, — ) =35.3023-sec ' -cm 



kef&ec6cl | Tavg , — ) = 0.7922* sec * 
\ "c^sc 2 / 

keflfrec6c7(T.z) 



•cm 



( 



I 



RatekC7H 1 4( 2- kptpdC3( T) , T)- Lcat / kgpbI6c7( T, z) 



RatekCH4(kpd(Tavg^^^j,Tavg^^j-Lcat=0.8I04'scc * -cm 

kgpb!6c7|Tavg^^.^Y5j = I0.1756-sec~' cm 



I 



keffrec6c7^Tavg^.^j =7.4067-sec ' "cm RatekC7H14(2-kptpdC3(Tavg^),Tavg^).Lcat =27.2196-sec" 



10 



kef&ec6cl (TrecZjj,zfec6jJ 



W — 

yj RatekCH4 (kpd (Treczj J , Trecz^ • Lcat q. 1 
|j kgpbl6cl(Trecz.^.zrrc6j 

Us 

a 



0.01 



ry 

j T=5 



no 



0.05 



Case specific ncase = 5 



10 



RatekCTH 1 4 i 2 kptpdC3 (Treczj^ , "^^^^iz) ' 
kgpbI6c7(Titcz.^, zrec6 . j 




0.1 
zrcc6. 



0.01 



0.15 



0.2 




Determine conversion of fuels with transport-limited finite catalytic combustion rates: 



Tree in Trecout 
Trecz. Trecin --""^ -"-^'^•iz 

iz ncasc „^ 



Case specific ncase = 5 



Trecz = 630.3722" K Trecout = 630.3722' K 

nz ncasc 



y2kc7hM- ^yc7hl40 



y2kch4Q "ych40 



yzkc7hl4.^^ J - [ 1 - keffrec6c7(Trecz.^,zrec6.J~ l.l- ^ 



£cell6 '^^^^'^^5 



, yzkc7iil4. 
dceli6 Urec Trecz. ^ 



iz/ 



yzkch4. =|l - keflTec6cl (Trecz. ,zrec6. 1.1 

izM I \ iz' iz/ 



ccelI6 



Trecin, 



, yzkch4. 

dceli6 Urec Trecz. / *^ 

ncase iz / 



Case specific 

ncase = 5 

Case specific 



5 ■ 
yj 



yzkc7hl4 =1.802' 10 



-5 



y2kch4^ = 0.0034 



yzkcThU-^ 
yzkcThM^ 0 1 

yzcThUQ 
yzkcM. 

~ 0.01 

yzkcb4f. 



no 



yzkc7h 14^=0.0005 
yzkch4jj= 0.005 



yzkc7hl4 yzkch4 

- 3.603- 10 ^ — — H 0.6838 



yzkc7hl4 



0.05 



O.I 
zrcc6. 



0.15 



0.2 



yzkch4Q 



Case specific ncase = 5 

Conversion of C7H14 under case 5 (full 
power) is 97.4% with estimated catalyfic 
rates and transport limttalions [vs. 99.6% 
for the maximum BL transport*limited 
conversion]. Conversion of CH4 is very 
poor (32%). CO conversion is prot>ably 
very high given higher catalytic and 
transport rates than with C7H14. 

This result shows that the specific rates 
of catalytic oxidation of all added 
components needs to be measured as 
applied to a recuperator. The use of a 
recuperator for emission control would 
appear to be greatly enhanced with 
higher exhaust (recuperator) 
temperatures. 



tn the following section, the conversion of C7H14, CH4 and CO Is considered for various honeycomb monliths 
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The. , Trecin. 

jcase t I jcasc t I 



pmixhc. . , pmixrecf Trecin. , ) 
jcasc ! I \ • jcase t I / 



Ghc2 



'jcase r 1 



mf 

Shcecell2 



Ghc4 



mf 

jcasc » I 



j^*^^* ShcecelW 



Ghc6. ^ , 

jcase -t- ! 



mf. 
jcasc 



Shc-ecell6 



Uhc2 = 



I 0 
3.058 
3.9753 
4.801 
5.8101 



m sec 



10.0353 



Uhc4 = 



0 

3.2157 
4.1804 
5.0487 
6.1098 
10.553 



m'sec 



Uhc6 = 



0 

3.1706 
4.1218 
4.9779 
6.0242 
10.4051 



* m* sec 



Evaluate Reynolds numbers for use in calculating the momentum, heat and mass transfer functions: 



S^^^^^^case., 



Q^'Sase-H -^^"^ 



Rehc4 



jcasc -I- 1 



Rehc6 



jcsse-f- 1 



Q^^^^jcase-Hl'^^^"^ 

^^'"*'^T(Thc.^^,j 



n 



tij 



ryRehc2 = 



0 

43.776 

56,9088 

68.7283 

83.1744 

167.4548 



Rehc4 = 



0 

31.7427 

41.2655 

49.836 

60.3111 

121.4242 



Rehc6 = 



0 

25.7355 
33.4562 
40.4048 
48.8975 
98.4453 



14.227 I 



dcell2 \ 



cfrichc2(j,z) . i I f 0.046263- f Rehc2. 

^ Rehc2. ! \ J z / 



0.45363 



cfrichc4a,z) -^^A I ^ 0.046263- (Rehc4.--—'^- 1 
^ Rehc4. ! \ J z / 

J - 



0.45363 



Average coefficients of friction 
for smooth square channels 



14.227 



dcell6 



, ;0.45363 



cfTichc6(j,z) -- - ■ • I r 0.046263-1 Rehc6. --- - ■ 
Rehc6. ' \ J z / ! 



CAPSTNE6.MCD l^HHH 
Calculate monolrth pressure drop in recuperator with entrance effect: 



='Uhc2. 

APhc2. . =4 cfrichc2(jcase i l.Lhc) - - ■ ^^^-t pmixhc . , 



Uhc4. 



I ^ . . Lhc 



APhc4. . =4 cfrichc4(jcasef- I.Lhc) J^?^± pmixhc. ^'"^ 



APhc6j^^i : = 4-cfrichc6(jcase^- l,Lhc)-^ pmixha^^^j- 



dcell6 





0 




10 




0 




0 




0.0013 




1 0.0028 




0.004 




0.0035 


APhc2 _ 


0.0017 


APhc4 


j 0.0036 


APhc6 _ 


0.0053 


3.5 APhc6 _ 


0.0046 


Phc 


0.002 


Phc 


i 0.0044 


Phc 


0.0064 


4 Phc 


0.0056 




0.0024 




1 0.0053 




0.0077 




0.0068 




. 0.0041 . 




[ 0.0088 . 




. 0.0127 . 




.0.0111 



APhcmax P = 



0 

0.003 
0.004 
0.005 
0.006 
0.011 



01 

a 

O 

03 

nJ 

ly 



Chose 600 cells per in2 honeycombs but 
reduce length to 3.5in via estimated AP 



Lhc = 3.5 in 



Use 600 cells/!n2 In honeycomb monolith for calculation of CH4, C7H14, and CO conversions 



Find mass transfer rates via Nusselt formulas: 



Niimavghc6cl .^^^ ^ i - 2.976 



I + 



7.6Lhc 



Rehc6. ^.ScITfThc. ^,Vdcell6l 

\ jcase f I V jcase 4- 1 y / 



Nuinavghc6cl = 



0 

2.976 

2.976 

2.9761 

2.9761 

2.9767 



. , ^ . / 7.6 Lhc 

Numavghc6co . . -2i976-! I ^ • : .77 

^ J«^«^^' lRehc6. -SccoTfThc. ,;.dcell6 

\ jcase t- 1 \ jcase 7 I ' 



Numavghc6co = 



0 

2.976 

2.976 

2.9761 

2.9761 

2.9769 



CAPSTNE6.MCD 



Numavghc6c7.^^^, =2.976 



I 

7.6'Lhc \ ^1 

Rehc6. ' * " •Sc7T(Thc" ^ , )-dcell6 

jcase -*- 1 \ jcasc -hi/ /J 



Nuniavghc6c7 



= i 



: 0 

! 2.9765 
' 2.977 
! 2.9777 
; 2.9788 
; 2.994 



Once again the entrance effects within monoliths are generally smalt 



Nowdetenmine the maximum possible CH4, C7H14. and CO conversions in a 600 cells per in2 monolith 

ycMO : = 0.005 yc7h 1 40 = 0.0005 ycoO = 0.005 



ych4hc6^^^^j — ych40-exp 



■ 4 • Lhc •Numavghc6c 1 



jcase + 1 



Rehc6. .SclTfThc. .Vdcell6 



jcasc -h I 



jcase -I- 1 /' 



- 4-Lhc-Numavghc6c 1 .^^ ^ , 

bj Rehc6. ^ -SclT/Thc ^,) dce!l6 
lease -t-l \ jcasc -1-1/ 





-64.0074 


Qi 


-49.2367 




-40.7696 




-33.689 


H 


- 16.7076 



/ 



4 • Lhc • N um avghc6co 



jcase -t- 1 



RJ ycohc6. ycoO' exp . , , 

^j-^ jcasci-i Rehc6 ^ -SccoT Thc. ^.)*dcelI6 

5 \ jcasc + I \ jcasc -H i / 



ycohc6 _ 
ycoO 



0 
0 
0 

j 4,7906- 10 
: 2.445- lO" 



ri4 



ych4hc6 _ 
ych40 



0 
0 
0 
0 

2.339- 10 



-15 



5.5463-10 



4 • Lhc -Num avghc6co 



jcase + \ 



Rehc6. 



jcase + 1 



SccoT.'Thc. . i-dcell6 
V jcase +- \ f 



58.27 



44.8234 



37.1153 



-30.6695 
- 15.2241 
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yc7hl4hc6. 



jcasc t I 



-yc7hI40exp 



yc7hl4hc6 
yc7hl40 



8.3588- iO 
3.5856- lO" 
3.0905-10" 

j 0.0002 

L 0.0141 



3 .5 • LhcNumavghc6c7 



j casc T 1 



jcase t I / 



dcelI6 



3.5LhcNumavghc6c7 



Rchc6. , •Sc7T(Thc. , /)-dcelI6 
jcasc i I \ jcasc t I / 



16.2974 



12.5386 



- 10.3846 



8.5841 



4.2631 



Only for C7H14 and only in case 5 (maximum power) does the BL transport-only 
calculated conversion (98.6%) fell below 99.9% with a 600 cel!&/in2 monolith With 
SO ppm C7H14 in the exhaust in this case, the exit concentration would be 0.7 ppm. 



^^fc Ift ASft A ft t ft ItH Ittit ItiMt ft ftft^ftft ftft ftft ft ftft ft f 



For the final analysis, consider the finite catalytic conversion of C7H1 4 and CH4 in a monolith: mease = 5 



Rehc6. ^ .Sc7T(Thc. ^,).dcell6 

jcasc -h 1 \ jcasc -f-l/ 



"43 zzcharhc6c7. ^ , 

,fl jcasc i-l 

hi 

m,^, , l^l^,.^.,SclT(T,.cl,^„).d»ll 6 

51 ^ ' 2 



Lhc 



zzcharhc6c7 



= 0.8137 



Lhc 



zzcharhc6cl 



= 2.8064 



,nz=400 



Cfl Lhc 

=3 ; zhc6. - iz* — -H dwall6 zzhc6c7 
« nz 

=•1 i 



Lhc = 8.89 'cm dcelI6 = 0.0961 -cm 
zhc6. 



Numono ^zzhc6c7j^ 
Nuinono(zzhc6cl|2) 



zhc6. 



zzcharhc6c7 



zzhc6cl. =- 



zzcharhc6c 1 



2zcharhc6c7 = 



mease 




0 

2.8414 
3.6938 
4.461 
5.3986 
[ 10.926 



"cm 



CAPSTNE6.MCD ^^^^^^^^^B 

I Case specific mease =5 

Dc7hl4T(T)Numono ^ 



Izzcharhc6c7 ' 

kgpblhc6c7(T,2) ^ ^ ^^^^ kgpblhc6c7(Thc .Lhc) = I3.6172'sec * "cm 

dceII6 ^ 

Dch4T(T)Numono( ^ 



I2zcharhc6cl j _ 

kgpblhc6cI(T,z) ^ ^ kgpblhc6cl(Thc ,Lhc) =46.9581-sec ' 'cm 

dcel!6 ^ ' 



kefflic6cl(T,z) 



I RatekCH4( kpd( T) . T)-Lcat / kgpblhc6c 1 ( T, z) 



Lhc \ 
^' 2 j 



kgpblhc6cl(Thc^^^.^^| =46.9589-sec ' 'cm 



_ keffhc6cl(Thc .—]= 1.8693 sec * cm 
^keffhc6c7(T.z) ^ 

j ; s 



RatekCH4(kpd(Thc ),Thc )-Lcat = 1.9467 "sec ' -cm 

\ V mease/ ' mease/ 



kgpbIhc6c7|Thc^,^j =13.6468-sec"' 'cm 



RatekC7H14{2-kptpdC3(T),T).Lcat/ kgpblhc6c7(T,z) 
keflBic6c7(Thc . — ] = 12.2537-sec'^ 'cm RatekC7H14(2-kptpdC3(Thc^J.Thc^^)-Lcat = 120.037rsec'' 'cm 



ffi 

m 
o 



10 



_R«ekC7Hi4 (2.kptpdC3 (Thc„^) .Tta^^) Let 
lti!pbl6c7(Thc^.ihc6i^. 



0.1 



















































































\ 
























F p 





0.02 0.04 0.06 0.08 

zhc6- 



0.1 
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Determine conversion of fuels with transport-limrted finite catalytic combustion rates: 



27 



Case specific mease = 5 



y2khc6c7hl4 yc7hl40 



yzkhc6ch4 - ych40 



yzkhc6c7hl4 ( I - keffhc6c7(Thc ,zhc6. ) ^^-^^ |.yzlchc6c7hl4. 

1 ^ nz dceII6 Uhc6 

^ mease/ 



Thc = 



i 0 

: 1005.3722 
I 1005.3722 
1005.3722 
1005.3722 
916.4833 



yzkhc6ch4 - ( I - keffhc6cl (Thc^ zhc6. ) -^^ -^^ 1 .yzkhc6ch4. 

'^^^ 1 ^ '^'^ nz dce!l6 Uhc6 ' 



Case specific 



mease =5 



yzkhc6c7hl4^ = 9.964- 10 



yzkhc6c7hl 4^=0.0005 



uJ 



y2khc6ch4^== 0.0028 



yzkhc6ch4^ = 0.005 



yzkhc6c7hl4 



yzkhc6c7hl4Q 



"^ = 1.993-10 ^ 



yzlchc6c7hl4j^ 
i^khc6c7hUQ 

y2khc6ch4. 



yzkhc6ch4Q 



no 




0.01 



yzkhc6ch4^ 
yzkhc6ch4 



= 0.5693 



Case specific mease = 5 

Conversion of C7H14 under case 5 (full 
power) with 600 cells/in2 x 3.5-ln length 
and extrapolated catalytic rates (via 
C3H8 at 400C) is 98.0% [vs. 99.6% via 
averaged integrated formulas, 99.3% 
with temperature dependent transport, 
and 98.6% with 400 cells/in2 and 4-in 
length. Conversion of CH4 is limited l>y 
slow Idnetics and CO is expected to t>e 
very high (>99.99%). 



zhc6. 
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Gas Properties and Constants 



*♦*♦**♦♦***♦**♦*♦*♦♦•*♦«************♦*****•■»***♦♦♦♦*♦*******■*♦♦*****#*** 



cm = 0.01m J^kg-m^-sec*^ 



-I 



W=Jsec 
atniE 101300- Pa 



Pref ^ l atm 
HmHO.OOOOOlm 



R:= 8.3143-:=^ 



Rc : = 8.3I43 J mole ' K 



-1 



n Gas 

1 CH4 

2 02 

3 N2 
4C02 
5H20 



MoIWl 



a(LJ) 



e(LJ) 



Define transport integrals: 
Constants: 

a^Hj =0.52487 



MWq := 16.044-gm Gq : = 3.785- 10' '^ m := 148.6 K 



MW, -3l!999 gin Cj : = 3.464- lO '^^ m e, = I06.7-K 



MWj : = 28.018-gni :^ 3.798- lO" *^-in -71.4-K 



MW3 : = 44.0l0 gm : = 3.94M0"'°-m S3 = 195.2'K 



MW4 := 18.015-gni = 2,641-10- '** m =809.1-K 



anbQ -0.14874 aDa^ = 1.06036 



CT^ib, =0.7732 CTDaj =0.193 



aDbp =0.1561 



oDb. =0.47635 



a^aj =2.16178 



Functions: 



oiib^ : = 2.43787 aDa^ : = 1 .03587 
0033= 1.76474 



Qy{T) = o^aQ-exp[ - (a^bQ-ln(T)) ] t ^ Cfiaj^ exp^- aubj^ T) 

1^=1 



aDbj- 1.52996 



aDbj =3.89411 



Qrt3) -1.0394 



OD(T) : ^ aDa(j-exp[ (aDbQ ln(T)) j -h ]^ aDajj^-exp^- aDbj^-T) 

iD= 1 



fiD(3) =0.95002 
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Set Inlet Conditions: 

P ;= l.O atm ych40 =0.035 

Define gas composition w/r methane using moist (3%) air; 



y(x)- 



( I - ych40) 0.2 - 2-ych40 -i- 2-x 
( 1 - ych40)-0.77 
ych40 - X 
( 1 - ych40) 0.03 + j2 ych40 - 2-x 



y(0.02) = 



0.02 

0.163 

0.743 

0.015 

0.059 



: a 

! . I 

i 



m 



m 



Define viscosity function 



2.6693- 10 



fi(T) = 



26 (MWo>T)°-^ io5.cn, /gm^o-^ 
'/ x2« (T\' sec ■Ik) 



2.6693- 10-^- 



10^ 



2.6693- 10 



2.6693 10 



26 



10^ 



sec 



cm] /gm 



cm I ^gm 

KT 



0.5 



0.5 



26/10^-cm\ /gm^^5 



l.ll6(MW^Ty 



0.5 



(3.115.10-«-cm)'.nji(^] 



2.6693 10^- 



10- 
sec 



'cm \ / gm 



0.5 



M(300-K) = 



1.103* 10 

2.06-10"^ 

1.77-10"^ 

1.5I9*10" 

1.068-10" 



•kg«m *»sec * 
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Hmix(x.T):= ^ 
n = 0 



y(x)„(n(T))„ 



4 

m = 0 


, /^'(T)„\°•'/Mw„^«•" 


2' 






0.5 



y(0.02) 



fiinix(0.02.300 K) = 1.759^10 ^•kg*m *-sec ' 

Define thermal conductivity functions 
igas =0,1 ..4 



0.02 

0.163 

0.743 

0.015 

0.059 



*********** 



1.869 




' 0.08727 


- .3273 




0.09966 


03919 


bX; = 


0.09816 


-7.215 




0.08015 


7.341 




-0.01013 



1.17910* 
-(3.743-10-^) 
-(5.067- 10-') 
5.477- 10-^ 
1.80M0'* 



(3.614- 10**) 
9.732- 10-' 

1.504- 10 * 
(1.053-10-*) 

(9.1-10-*) 



0.00001 



w 

ctn-K 



W/cm/K 



Xmix(x,T) := ^ 
i = 0 



y(xVX(i.T) 



4 

j = o 




2 




r / MW,\i 


0,5 



CH4 

02 
N2 

C02 

H20 



X< igas, 1000- K) 


0,16716-kg-m-sec'^ 




0.07163-kgin-sec'^ 




0.06292kg-m-sec"^ 




0.06788 kg-ni-sec^ 




0.0863l-kg-m-sec^ 





Xmix(.0, lOOO K) = 6.64- 1 0"** ~ 

cm-K 



********************««**««********«**««*«*■* 



# 



**************************************************************************************************************** 



Define thermal conductivity of alumina and SS321 



Altc 



r 83.71 
I 224.898 
255.842 
132.727 
25.9695 



The parameter set Altc_n was determined from published values of the 
thermal conductivity, Ref: HCP 

AAI203(x) is the thermal conductivity function (of x with unit K) while 
XI is its derivative w/r to temperature. 



XAI203(x):. S ToOOlC STK ^*('^> = E TOOOTR-"' TOOa 

n=0 ^ ' n=l ^ 



n-I 



W 
nTK 



S32ltc : 



r 4.69312 
4.9812 
0.33357 

0 

0 



The parameter set S321tc_n was determined from published values of the 
thermal conductivity. Ref: Handbook of Metals 8th Ed. Vol. 1 p. 423 



;iss32l(x) := S321tc„ 
n = 0 



/ X \" W 



01 



W W 
Xss32 1 ( 1 OOOK) = 9.34 1 *^=. >AI203( 1 OOOK) = 7.897 ^ 



i =30,31.. 200 



Ta::=IOiK 



W 



\1{ lOOOKVlOOOK =-7.517 



m 
□ 



40 



30 



XAl203^aj) 20 
Xl^aj^lOOaK 
^321 (Ta.) ,0 



-I 
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X M500 K xO root(XI{x).x) 



xO = 1.588MO^K 



Xss321(300 K) =6.157- 



W 



O.OOS in =7.62-10 ^-cm 



XA12O3(300-K) =0.359' 



W 
cmK 



XAI203( 200- K) =0.479" 



W 
cmK 



XAI2O3(800 K) =0.102- 



W 
cm-K 



^I203(x0) =0.054- 



W 



Minimum Ihemial 



cm-K conductivity 



Detemnine themnal conductivity of porous wall: 

Ref: RHG - Handbook Chpt 6 for effective themial conductivity including 
convection and radiative terms 



erad:-0.6 e:=0.42 



dp : = 0.00005 cm 



dp =5-10~^Tn 



+ 0.229- 10*— ^e-dperadT^le + J— ^ 



0.229- 10*— .j—je-dp-erad-T^ 
XA1203(T) 



^ ^ Radiation component 

Xwc(0.01,1000 K) =o.473--^-j^ Xwcnr(0.0l . lOOO K) =0.473-^ negligible with such small 

particles In the washcoat 



w w 

Xmix(O.OMOOOK) =0.067- MI2O3(700 K) = 12.354---. 

m-K. m*K. 



Xss321(IOOOK) =9.341- 



W 
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XAI2O3(1000 K) 
Xmix(O.OMOOOK) 



117.74 



Xwc(0,01,1000K) 
Xmix(O.Ol.lOOOK) 



The thermal conductivity of the wash 
= 7.046 coat is about 7 times greater at 1 0OOK 
than the gas phase conductivity 



Xwc(O.OMOOOK) 
Ais321(1000K) 



= 0.051 



Xwc(O.0l, lOOO K) Xss321(IOOO K) 



0.001 cm 



0.003- in 



Xss32I(1000K) 
0.003 in 



= 1.386 



The washcoat increases 
the thenmal conductance of 
the wall by about 39% 



Xss321( 1000 K) 
Xmix(O.Ol.lOOOK) 



139.275 



For a Nusselt no. of 2.976 (dcellS = 0.0961 *cm) for developed flow in the recuperator (full power) the 
mean heat transfer boundaiSy layer thickness is much greater than either the SS321 wall thickness or 
the washcoat thickness. 



Sgprecup : 



0.0961 cm 
2.976 



6gprecup = 0.032 -cm 



hrecavg(T) : = 



/ ^ \ 


1 


Amix(0.01,T) 


^ Xss321(1000 K) 



Sgprecup / 



0.003-in 



Overall local heat transfer coefficient 
for the recuperator (assuming gas phase 
conductivities are about the same 



hrecavgwc(T) 



Effect of wash coat on heat 



1 




/ ^ \ 


1 transfer doefficient 


Ajnix(0.01,T) 


Xmix(0.01,T) 


/Xss32l(1000 K)\ ' A.wc(0.01.T) 


20.001cm 
Sgprecup-- 2^^^ 




\ figprccup 1 


\ 0.003in / 0.001-cm 



hrecavg( 1 000' K)= 0.0 1038 



W 



cm^K 



hrecavgwc(1000-K) =0.01046- 



W 



The wash coat makes about 0.1% increase in heat transfer 
coefficient because of the slightly decreased cell size and greater 
thermal conductivity of the washcoat vs. gas phase. 



